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20 BACKGROUND 

Field of the Invention 

The present invention relates to performing arithmetic operations on 
interval operands within a computer system. More specifically, the present 
25 invention relates to a method and an apparatus for finding the roots of a nonlinear 
equation of one variable using the interval version of Newton's method. 

Related Art 

Rapid advances in computing technology make it possible to perform 
30 trillions of computational operations each second. This tremendous 
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computational speed makes it practical to perform computationally intensive tasks 
as diverse as predicting the weather and optimizing the design of an aircraft 
engine. Such computational tasks are typically performed using machine- 
representable floating-point numbers to approximate values of real numbers. (For 

5 example, see the Institute of Electrical and Electronics Engineers (IEEE) standard 
754 for binary floating-point numbers.) 

In spite of their limitations, floating-point numbers are generally used to 
perform most computational tasks. 

One limitation is that machine-representable floating-point numbers have a 

10 fixed-size word length, which limits their accuracy. Note that a floating-point 
number is typically encoded using a 32, 64 or 128-bit binary number, which 
means that there are only 2 32 5 2 64 or 2 128 possible symbols that can be used to 
specify a floating-point number. Hence, most real number values can only be 
approximated with a corresponding floating-point number. This creates 

1 5 estimation errors that can be magnified through even a few computations, thereby 
adversely affecting the accuracy of a computation. 

A related limitation is that floating-point numbers contain no information 
about their accuracy. Most measured data values include some amount of error 
that arises from the measurement process itself. This error can often be quantified 

20 as an accuracy parameter, which can subsequently be used to determine the 

accuracy of a computation. However, floating-point numbers are not designed to 
keep track of accuracy information, whether from input data measurement errors 
or machine rounding errors. Hence, it is not possible to determine the accuracy of 
a computation by merely examining the floating-point number that results from 

25 the computation. 

Interval arithmetic has been developed to solve the above-described 
problems. Interval arithmetic represents numbers as intervals specified by a first 
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(left) endpoint and a second (right) endpoint For example, the interval [a, b] 9 
where a < b, is a closed, bounded subset of the real numbers, R, which includes a 
and b as well as all real numbers between a and b. Arithmetic operations on 
interval operands (interval arithmetic) are defined so that interval results always 

5 contain the entire set of possible values. The result is a mathematical system for 
rigorously bounding numerical errors from all sources, including measurement 
data errors, machine rounding errors and their interactions. (Note that the first 
endpoint normally contains the "infimum", which is the largest number that is less 
than or equal to each of a given set of real numbers. Similarly, the second 

1 0 endpoint normally contains the "supremum", which is the smallest number that is 
greater than or equal to each of the given set of real numbers.) 

One commonly performed computational operation is to find the roots of a 
nonlinear equation using Newton's method. The interval version of Newton's 
method works in the following manner. From the mean value theorem, 

15 

f(x)-f(x*) = (x-x*)fW> 

where £ is some generally unknown point between x and x* If x* is a zero of/ 
then f(x*) = 0 and, from the previous equation, 

20 

x*=x-f(x)/fW- 

Let X be an interval containing both x and x* Since f is between x and x* it 
follows that £ eX. Moreover, from basic properties of interval analysis it follows 
25 that/Y# efW Hence, x* eN(x,X) where 

N(xJ)=x-f(x)/f'(X). 

3 

Attorney Docket No. SUN-P6282-SPL Inventors: Waister, et al. 

ARPHASUN MICROSYSTEMS\SUN-P6282-SPL\SUN-P6282-SPL APPLICATION3.DOC 



Temporarily assume 0 $f*(X) so that N(x, X) is a finite interval Since any zero of 
/in X is also in N(x,X) 9 the zero is in the intersection X n N(x y X). Using this fact, 
we define an algorithm for finding zero x* LetXo be an interval containing jc* 
5 For n = 0, 1, 2, define 

X n = m(X tl ) 

N(x n) X n )=x n -f(x n )/f'(X n ) 
X n+ i =X n nN(x m Xr), 

10 

wherein m(X) is the midpoint of the interval X. We call x n the point of expansion 
for the Newton method. It is not necessary to choose x n to be the midpoint of X n . 
The only requirement is that x n e X n to assure that x * e N(x n ,X r} ). However, it is 
convenient and efficient to choose x„ = m(X n ). 
1 5 Roots of an interval equation can be intervals rather than points when the 

equation contains non-degenerate interval constants or parameters. Suppose the 
interval version of Newton's method to find the roots of a single nonlinear 
equation has not yet satisfied the user-specified convergence tolerances. Then it is 
difficult to distinguish between the following three situations: 

20 

a) the current interval is a tight enclosure of a single interval root; 

b) the current interval contains sufficiently distinct interval roots that they 
can be isolated with a reasonable amount of effort; and 

c) the current interval contains point and/or interval roots that are so close 
25 to being indistinguishable that the effort to isolate them is unreasonably 

large with the existing wordlength. 
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What is needed is a method and an apparatus for terminating the interval 
version of Newton's root finding method before iterations lose their practical 
value in isolating meaningfully distinct interval roots. 

5 SUMMARY 

One embodiment of the present invention provides a system for finding 
zeros of a function,/ within an interval, X 9 using the interval version of Newton's 
method. The system operates by receiving a representation of the interval X This 
representation including a first floating-point number, a, representing the left 

1 0 endpoint of X, and a second floating-point number, &, representing the right 

endpoint of X. Next, the system performs an interval Newton step onX, wherein 
the point of expansion is the midpoint, x, of the interval X. Note that performing 
the interval Newton step involves evaluating /fc) to produce an interval result 
f J (x). Iff^x) contains zero, the system evaluates/^ to produce an interval result 

15 f l ( a )- It also evaluates/^ to produce an interval result/ 7 ^. The system then 
evaluates a termination condition for the processing of the current interval X, 
wherein the termination condition is TRUE if a number of sub-conditions are 
satisfied, including iff 1 (a) contains zero and if f J (b) contains zero. If the 
termination condition is TRUE, the system terminates the processing of the 

20 current interval X 3 and records X as a final bound. 

In one embodiment of the present invention, if f l (a) does not contain zero, 
evaluating f(a) additionally involves performing an interval Newton step wherein 
the point of expansion is a. 

In one embodiment of the present invention, if f ! (b) does not contain zero, 

25 evaluating/^ additionally involves performing an interval Newton step wherein 
the point of expansion is b. 
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In one embodiment of the present invention, iff(x) contains zero and if 
/ '(X) contains zero, the termination condition for processing the current interval X 
is TRUE iff 1 (a) contains zevoj 1 ^) contains zero, f ! (xj) contains zero and if 
f l (x 2 ) contains zero. Note that xj is the midpoint between a and x; and x 2 is the 
5 midpoint between x and b. 

In one embodiment of the present invention, if f(x) contains zero and if 
f'(X) contains zero, and if either f(x } ) or f(x 2 ) does not contain zero, the system 
additionally splits the interval Xin half and applies the interval Newton method to 
each half separately. 

10 In one embodiment of the present invention, in the case where f J (x) does 

not contain zero and if f'(X) contains zero, if the width of the interval X divided 
by the magnitude of the interval X is less than a first threshold value, and the 
magnitude of f(X) is less than a second threshold value, the method further 
comprises terminating the processing of the current interval^ and recording X as 
1 5 a final bound. 

In one embodiment of the present invention, if a given Newton step does 
not reduce the width of an interval by at least half, the system additionally splits 
the interval in half and applies the interval Newton method to each half separately. 
In one embodiment of the present invention, if an interval Newton step 
20 results in two intervals, the system additionally applies the interval Newton 
method to each of the two intervals separately. 

In one embodiment of the present invention, if the result of an interval 
Newton step is the empty interval, the system returns to process another interval. 
In one embodiment of the present invention, iff 1 (a) contains zero and if 
25 f J (b) contains zero, the system returns to process another interval. 

In one embodiment of the present invention, if f(a) does not contain zero 
or if f l (b) does not contain zero, the system performs an interval Newton step 
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wherein the point of expansion is the midpoint of the interval X If the result of 
this interval Newton step is the empty interval, the system returns to process 
another interval. 

In one embodiment of the present invention, the system terminates the 
5 processing of the current interval X after a predetermined number of iterations. 

BRIEF DESCRIPTION OF THE FIGURES 

FIG. 1 illustrates a computer system in accordance with an embodiment of 
the present invention. 
10 FIG. 2 illustrates the process of compiling and using code for interval 

computations in accordance with an embodiment of the present invention. 

FIG. 3 illustrates an arithmetic unit for interval computations in 
accordance with an embodiment of the present invention. 

FIG. 4 is a flow chart illustrating the process of performing an interval 
1 5 computation in accordance with an embodiment of the present invention. 

FIG. 5 illustrates four different interval operations in accordance with an 
embodiment of the present invention. 

FIG. 6 illustrates the process of solving for zeros of a function using the 
interval Newton method in accordance with an embodiment of the present 
20 invention. 

FIG. 7 illustrates the interval Newton method for the case where the 
derivative of the function does not contain zero in accordance with an 
embodiment of the present invention. 

FIG. 8 illustrates the interval Newton method for the case where the 
25 derivative of the function contains zero in accordance with an embodiment of the 
present invention. 
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DETAILED DESCRIPTION 

The following description is presented to enable any person skilled in the 
art to make and use the invention, and is provided in the context of a particular 
5 application and its requirements. Various modifications to the disclosed 

embodiments will be readily apparent to those skilled in the art, and the general 
principles defined herein may be applied to other embodiments and applications 
without departing from the spirit and scope of the present invention. Thus, the 
present invention is not intended to be limited to the embodiments shown, but is 
10 to be accorded the widest scope consistent with the principles and features 
disclosed herein. 

The data structures and code described in this detailed description are 
typically stored on a computer readable storage medium, which may be any device 
or medium that can store code and/or data for use by a computer system. This 

1 5 includes, but is not limited to, magnetic and optical storage devices such as disk 
drives, magnetic tape, CDs (compact discs) and DVDs (digital versatile discs or 
digital video discs), and computer instruction signals embodied in a transmission 
medium (with or without a carrier wave upon which the signals are modulated). 
For example, the transmission medium may include a communications network, 

20 such as the Internet. 

Computer System 

FIG. 1 illustrates a computer system 100 in accordance with an 
embodiment of the present invention. As illustrated in FIG. 1, computer system 
25 100 includes processor 102, which is coupled to a memory 1 12 and a to peripheral 
bus 1 10 through bridge 106. Bridge 106 can generally include any type of 
circuitry for coupling components of computer system 100 together. 
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Processor 102 can include any type of processor, including, but not limited 
to, a microprocessor, a mainframe computer, a digital signal processor, a personal 
organizer, a device controller and a computational engine within an appliance. 
Processor 102 includes an arithmetic unit 104, which is capable of performing 
5 computational operations using floating-point numbers. 

Processor 102 communicates with storage device 108 through bridge 106 
and peripheral bus 110. Storage device 108 can include any type of non- volatile 
storage device that can be coupled to a computer system. This includes, but is not 
limited to, magnetic, optical, and magneto-optical storage devices, as well as 
10 storage devices based on flash memory and/or battery-backed up memory. 

Processor 102 communicates with memory 1 12 through bridge 106. 
Memory 1 12 can include any type of memory that can store code and data for 
execution by processor 102. As illustrated in FIG. 1, memory 1 12 contains 
computational code for intervals 1 14. Computational code 114 contains 
1 5 instructions for the interval operations to be performed on individual operands, or 
interval values 115, which are also stored within memory 112. This 
computational code 1 14 and these interval values 1 15 are described in more detail 
below with reference to FIGs. 2-5. 

Note that although the present invention is described in the context of 
20 computer system 100 illustrated in FIG. 1, the present invention can generally 

operate on any type of computing device that can perform computations involving 
floating-point numbers. Hence, the present invention is not limited to the 
computer system 100 illustrated in FIG. 1. 



25 Compiling and Using Interval Code 

FIG. 2 illustrates the process of compiling and using code for interval 
computations in accordance with an embodiment of the present invention. The 
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system starts with source code 202, which specifies a number of computational 
operations involving intervals. Source code 202 passes through compiler 204, 
which converts source code 202 into executable code form 206 for interval 
computations. Processor 102 retrieves executable code 206 and uses it to control 
5 the operation of arithmetic unit 1 04. 

Processor 102 also retrieves interval values 115 from memory 1 12 and 
passes these interval values 1 15 through arithmetic unit 104 to produce results 
212. Results 212 can also include interval values. 

Note that the term "compilation" as used in this specification is to be 
1 0 construed broadly to include pre-compilation and just-in-time compilation, as well 
as use of an interpreter that interprets instructions at run-time. Hence, the term 
"compiler" as used in the specification and the claims refers to pre-compilers, 
just-in-time compilers and interpreters. 



15 Arithmetic Unit for Intervals 

FIG. 3 illustrates arithmetic unit 104 for interval computations in more 
detail accordance with an embodiment of the present invention. Details regarding 
the construction of such an arithmetic unit are well known in the art. For 
example, see U.S. Patent Nos. 5,687,106 and 6,044,454, which are hereby 
20 incorporated by reference in order to provide details on the construction of such 
an arithmetic unit. Arithmetic unit 104 receives intervals 302 and 312 as inputs 
and produces interval 322 as an output. 

In the embodiment illustrated in FIG. 3, interval 302 includes a first 
floating-point number 304 representing a first endpoint of interval 302, and a 
25 second floating-point number 306 representing a second endpoint of interval 302. 
Similarly, interval 312 includes a first floating-point number 314 representing a 
first endpoint of interval 312, and a second floating-point number 316 
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representing a second endpoint of interval 312. Also, the resulting interval 322 
includes a first floating-point number 324 representing a first endpoint of interval 
322, and a second floating-point number 326 representing a second endpoint of 
interval 322. 

5 Note that arithmetic unit 1 04 includes circuitry for performing the interval 

operations that are outlined in FIG. 5. This circuitry enables the interval 
operations to be performed efficiently. 

However, note that the present invention can also be applied to computing 
devices that do not include special-purpose hardware for performing interval 

10 operations. In such computing devices, compiler 204 converts interval operations 
into a executable code that can be executed using standard computational 
hardware that is not specially designed for interval operations. 

FIG. 4 is a flow chart illustrating the process of performing an interval 
computation in accordance with an embodiment of the present invention. The 

1 5 system starts by receiving a representation of an interval, such as first floating- 
point number 304 and second floating-point number 306 (step 402). Next, the 
system performs an arithmetic operation using the representation of the interval to 
produce a result (step 404). The possibilities for this arithmetic operation are 
described in more detail below with reference to FIG. 5. 

20 

Interval Operations 

FIG. 5 illustrates four different interval operations in accordance with an 
embodiment of the present invention. These interval operations operate on the 
intervals X and 7. The interval X includes two endpoints, 

x denotes the lower bound of X, and 
x denotes the upper bound of X, 

25 
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The interval X is a closed subset of the extended (including -°° and +«) real 
numbers i?* (see line 1 of FIG. 5). Similarly the interval 7 also has two endpoints 
and is a closed subset of the extended real numbers i?* (see line 2 of FIG. 5). 

Note that an interval is a point or degenerate interval if X= [x, x]. Also 
5 note that the left endpoint of an interior interval is always less than or equal to the 
right endpoint. The set of extended real numbers, R* is the set of real numbers, R 9 
extended with the two ideal points negative infinity and positive infinity: 

R* = Ru {-oo} u {+oc}. 

10 

In the equations that appear in FIG. 5, the up arrows and down arrows 
indicate the direction of rounding in the next and subsequent operations. Directed 
rounding (up or down) is applied if the result of a floating-point operation is not 
machine-representable. 
1 5 The addition operation X+ 7 adds the left endpoint of X to the left 

endpoint of 7 and rounds down to the nearest floating-point number to produce a 
resulting left endpoint, and adds the right endpoint of AT to the right endpoint of Y 
and rounds up to the nearest floating-point number to produce a resulting right 
endpoint. 

20 Similarly, the subtraction operation X - Y subtracts the right endpoint of Y 

from the left endpoint of X and rounds down to produce a resulting left endpoint, 
and subtracts the left endpoint of 7 from the right endpoint of X and rounds up to 
produce a resulting right endpoint. 

The multiplication operation selects the minimum value of four different 

25 terms (rounded down) to produce the resulting left endpoint. These terms are: the 
left endpoint of X multiplied by the left endpoint of 7; the left endpoint of X 
multiplied by the right endpoint of 7; the right endpoint of X multiplied by the left 
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endpoint of 7; and the right endpoint of X multiplied by the right endpoint of Y. 
This multiplication operation additionally selects the maximum of the same four 
terms (rounded up) to produce the resulting right endpoint. 

Similarly, the division operation selects the minimum of four different 

5 terms (rounded down) to produce the resulting left endpoint. These terms are: the 
left endpoint of X divided by the left endpoint of 7; the left endpoint of X divided 
by the right endpoint of 7; the right endpoint of X divided by the left endpoint of 
7; and the right endpoint of X divided by the right endpoint of 7. This division 
operation additionally selects the maximum of the same four terms (rounded up) 

10 to produce the resulting right endpoint. For the special case where the interval 7 
includes zero, X/Yis an exterior interval that is nevertheless contained in the 
interval i?*. 

Note that the result of any of these interval operations is the empty interval 
if either of the intervals, Xor 7, are the empty interval. Also note, that in one 
1 5 embodiment of the present invention, extended interval operations never cause 
undefined outcomes, which are referred to as "exceptions" in the IEEE 754 
standard. 

Interval Version of Newton's Method 

20 FIG. 6 illustrates the procfess of solving for zeros of a function using the 

interval version of Newton's method in accordance with an embodiment of the 
present invention. The system starts with an initial interval X 0 and stopping 
tolerances s x and s F . The system stops when interval bounds for all zeros of fin 
X 0 have been found. After termination, the bounds on any simple zero generally 

25 approximate the optimal bound. If the tolerances s x and s F are not chosen too 

small, then each multiple zero of /in X 0 is generally isolated within an interval of 
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relative width less than Also, we generally have \f(x)\ < s F for all points x in an 
interval bounding a zero. 

In the following description, the current interval is denoted by Xat each 
step although X changes from step to step. Also, x denotes the midpoint m(X) so x 
5 changes as well. Note that except for a few cases described below, the Newton 
step is generally defined by an expansion about x. 

The system starts by placing the initial interval, X#, in a list of intervals to 
be processed, L (step 602). 

If L is empty, the system stops. Otherwise, the system selects the interval 
1 0 X, from L that has been in L for the shortest amount of time, and the system 
deletes X from L (step 604). 

Next, the system determines if the derivative of/ which is represented as 
when evaluated over the interval X contains zero. In other words, is 0 ef'(X)l 
(step 606). If not, the system finds roots for the case where f'(X) does not contain 
1 5 zero (step 608). The operations involved in solving for this case are described in 
more detail with reference to FIG. 7 below. The system then returns to step 604 
process additional intervals if necessary. 

lff'(X) contains zero at step 606, the system determines whether f\x) 
contains zero (step 610). Note that/ 7 ^ denotes the interval resulting from an 
20 evaluation of/fx). If the conditions in both steps 606 and 610 are true, the system 
finds the roots of /for the case where f'(X) does not contain zero (step 612). The 
operations involved in solving for this case are described in more detail with 
reference to FIG. 8 below. The system then returns to step 604 to process 
additional intervals if necessary. 
25 K/ 7 fo) does not contain zero at step 610, the system performs tests 

involving the stopping tolerances s x and s F (step 614). If the width of the interval 
X (denoted as w(X)) divided by the magnitude of X (denoted as |X|)is less than s x> 
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and if the magnitude offlJC) (denoted as \f{X)\) is less than e F , the system records X 
as a final bound (step 616) and returns to step 604 to process additional intervals 
if necessary. Note that the width of the interval [a 9 b] (denoted as w(X)) is 
simply b - a. Also note that the magnitude of an interval X= [a, b] (denoted by 

5 \X\) is the maximum value of |x| for all x e X. Thus, \X\ = max(|a|,|6|). 

Furthermore, note that the stopping tolerances, Sx and 6f, are used to stop the 
system from performing further computations when the interval bound on a zero 
of /is sufficiently narrow. 

If at step 614 one or both of the tested quantities is greater than or equal to 

10 the corresponding stopping tolerance, the system applies the interval Newton step 
to the interval X (step 617). 

The system then examines the result of the Newton step (step 618). If the 
result is empty, the system returns to step 604 to process additional intervals, if 
necessary. If the result is a single interval, the system determines whether the 

1 5 Newton step reduced the interval X by at least half in size (step 620). If so, the 
system returns to step 606 to continue processing the current interval X. 
Otherwise, the system splits X into halves, puts one half into the list of intervals to 
be processed, L, and designates the other half as the current interval X (step 622). 
The system then returns to step 606 to continue processing the current interval X. 

20 If the result of the Newton step at step 618 is two intervals, the system puts 

one of the two intervals into the list of intervals to be processed, L 9 and designates 
the other interval as the current interval, X (step 624). The system then returns to 
step 606 to continue processing the current interval X. 

25 Case Where 0 ef'OQ 

FIG. 7 illustrates the interval Newton method for the case where the 
derivative of the function does not contain zero in accordance with an 

15 

Attorney Docket No. SUN-P6282-SPL Inventors: Walster, et al. 

ARPHASUN MICROSYSTEMS\SUN-P6282-SPL\SUN-P6282-SPL APPLICATION3.DOC 



embodiment of the present invention. This flow chart illustrates in more detail the 
operations carried out at step 608 of FIG. 6. 

Note that the relation 0 gf'QQ assures that there is no more than one point 
or interval zero of/in X Moreover, if there is a zero of /in AT, then it is a simple 
5 one. Also note that inclusion isotonicity assures that 0 gf'(X') for any interval X' 
contained inX. 

The system starts by initializing a counter variable, n, and two flags, F a 
and Ffo to zero (step 702). Note that the system sets F a = 1 when 0 ef'(a). 
Similarly, the system sets F b = 1 when 0 ef'(b). Also note that the system cycles 
1 0 through the operations listed below no more than four times. The integer n counts 
the cycles. 

After initializing the variables in step 702, the system tests F a to see if 
F a ==l (step 703). If F a == h the system proceeds to step 710 to test flag F h . On 
the other hand, if F a * 1 at step 703, the system evaluates f(a) to produce the 

15 resulting interval f(a) (step 704). Next, the system determines if f(a) contains 
zero (step 706). If so, the system sets F a =l (step 707) and proceeds to step 71 0 to 
test flag F b . If at step 706 f(a) does not contain zero, the system applies an 
interval Newton step wherein the point of expansion is "a" (step 708). If the 
result of this expansion is the empty interval, the system returns to process the 

20 next interval. 

Otherwise, the system tests F b to see if F b ==l (step 710). If F*==i 5 the 
system proceeds to step 7 1 8 to test both flags F a and F b . On the other hand, if 
F b #7 at step 710, the system evaluates f(b) to produce the resulting interval/ 7 ^ 
(step 712). Next, the system determines if f(b) contains zero (step 714). If so, 

25 the system sets F b =l (step 7 1 5) and proceeds to step 7 1 8 to test both flags F a and 
F b . If at step 714 f(b) does not contain zero, the system applies an interval 
Newton step wherein the point of expansion is "6" (step 716). If the result of this 
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expansion is the empty interval, the system returns to process the next interval. 
Let x** be the midpoint of X that satisfies the condition 0 ef(x) in step 610. x** 
is set in step 612. Note that if x** is no longer an element of X, the system uses 
the midpoint of the currents m(X), as the point of expansion. 
5 Otherwise, the system tests to see if both F a ==l and F b ==l (step 718). If 

both F a ==l and F*==i, the bounds, a and b, on the zero of / are as tight as they 
are likely to become. Hence, the system returns to process the next interval (step 
721). 

Otherwise, the system applies an interval Newton step wherein the point of 
10 expansion is m(X) (step 720). If the result of this expansion is the empty interval, 
the system returns to process the next interval. Next, the system increments n 
(step 722) and tests to see if n < 4 (step 724). If n > 4 } the system returns to 
process the next interval because the bounds on the zero are unlikely to be further 
tightened by additional processing (step 726). Otherwise, if n < 4, the system 
1 5 returns to step 703 to repeat the process. 

Case Where 0 ef VC) 

FIG. 8 illustrates the interval Newton method for the case where the 
derivative of the function contains zero and 0 ef(x) in accordance with an 
20 embodiment of the present invention. This flow chart illustrates in more detail the 
operations carried out at step 612 of FIG. 6. 

The system first evaluates f(a) to produce the resulting interval/ 7 ^ 
(step 802). Next, the system determines if f(a) contains zero (step 804). If not, 
the system applies an interval Newton step wherein the point of expansion is "a" 
25 (step 806) and returns to the main program. 

On the other hand, if f(a) contains zero, the system evaluates/^ to 
produce the resulting interval f(b) (step 808). Next, the system determines if 
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f(b) contains zero (step 810). If not, the system applies an interval Newton step 
wherein the point of expansion is "6" (step 812) and returns to the main program. 

On the other hand, if f(b) contains zero, the system next finds the centers, 
x; and x 2 , of the two halves of the interval X by computing x } = (3a + b)/4 and 

5 X2 = ( a + sb)/4 (step 814). The system then evaluates the function/at each of the 
centers, x } and x 2 , to produce the resulting intervals f(%i) mdf ! (x 2 ). Next, the 
system determines if both f(xj) contains zero and f(x 2 ) contains zero (step 816). 
If so, the computed (interval) value of / contains zero for each of the five equally 
spaced points, a, x h m(X) 9 x 2 and b. In this case, additional processing is unlikely 

1 0 to narrow the X further. Hence, the system accepts X as a final bound on a zero of 
/(step 822). 

On the other hand, if either f J (x } ) or f(x 2 ) does not contain zero, the 
system splits the interval X in half, and applies the interval Newton method to 
each half separately (step 820). 

1 5 Note that the system can also make use of higher precision interval 

arithmetic to help resolve the uncertainty when 0 e f(x) and 0 e f'(X). 

The foregoing descriptions of embodiments of the present invention have 
been presented for purposes of illustration and description only. They are not 
intended to be exhaustive or to limit the present invention to the forms disclosed. 

20 Accordingly, many modifications and variations will be apparent to practitioners 
skilled in the art. 

For example, in the computational process described above with reference 
to FIG. 6, it is possible to process intervals in the list, Z, in parallel using multiple 
concurrently executing processors. Furthermore, it is possible to subdivide the 
25 initial interval into multiple smaller sub-intervals so that the multiple sub-intervals 
can be processed in parallel 
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In another embodiment of the present invention, if the function,/ has 
already been evaluated at a right endpoint or a left endpoint of an interval, X, the 
system saves the result of this evaluation in the list, L 9 so that the function does 
not have to be reevaluated at the right endpoint or left endpoint of X. 

Hence, the above disclosure is not intended to limit the present invention. 
The scope of the present invention is defined by the appended claims. 



19 

Attorney Docket No. SUN-P6282-SPL Inventors: Walster, et al. 

ARPH:\SUN MICR0SYSTEMS\SUN-P6282-SPL\SUN-P6282-SPL APPLICATIONS DOC 



